// Magic Software, Inc.
// http://www.magic-software.com
// Copyright (c) 2000, All Rights Reserved
//
// Source code from Magic Software is supplied under the terms of a license
// agreement and may not be copied or disclosed except in accordance with the
// terms of that agreement.  The various license agreements may be found at
// the Magic Software web site.  This file is subject to the license
//
// FREE SOURCE CODE
// http://www.magic-software.com/License/free.pdf

#ifndef MGCMIDPOINT_H
#define MGCMIDPOINT_H

#include "MgcODE.h"


class MgcMidpoint : public MgcODE
{
public:
    virtual ~MgcMidpoint ();

    // for dx/dt = F(t,x)
    MgcMidpoint (int iDim, MgcReal fStep, Function* aoF);
    virtual void Update (MgcReal fTIn, MgcReal* afXIn, MgcReal& rfTOut,
        MgcReal* afXOut);

    // for dx/dt = A(x)
    MgcMidpoint (int iDim, MgcReal fStep, AutoFunction* aoA);
    virtual void Update (MgcReal* afXIn, MgcReal* afXOut);

    virtual void SetStepSize (MgcReal fStep);

protected:
    MgcReal m_fHalfStep;
    MgcReal* m_afXTemp;
};

#endif
